o 

O 
< 



Quantum electrodynamics and photon-assisted tunnelling in long Josephson junctions 

A.O. Sboychakov,^'^ Sergey Savel'ev,^''^ and Franco Nori^'^ 

' Frontier Research System, The Institute of Physical and Chemical Research (RIKEN), Wako-shi, Saitama, 351-0198, Japan 

''institute for Theoretical and Applied Electrodynamics Russian 

Academy of Sciences, Izhorskaya str. 13, Moscow 125412, Russia 

Department of Physics, Loughborough University, Loughborough LEll 3TU, United Kingdom 

^MCTP, CSCS, Department of Physics, University of Michigan, Ann Arbor, MI 48109, USA 

(Dated: April 14, 2008) 

We describe the interaction between an electromagnetic field and a long Josephson junction (JJ) 
driven by a dc current. We calculate the amplitudes of emission and absorption of light via the 
creation and annihilation of quantized Josephson plasma waves (JPWs). Both, the energies of JPW 
quanta and the amplitudes of light absorption and emission, strongly depend on the junction's 
length and can be tuned by an applied dc current. Moreover, photon-assisted macroscopic quantum 
tunnelling in long Josephson junctions show resonances when the frequency of the outside radiation 
coincides with the current-driven eigenfrequencies of the quantized JPWs. 
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I. INTRODUCTION 

The miniaturization of electronic devices allows the ob- 
servation of quantum effects which were impossible to 
measure in the past. Systems of Josephson junctions, 
characterized by high frequency (up to several THz), ex- 
hibit a crossover to the quantum regime at relatively high 
temperatures. Indeed, quantum oscillations and macro- 
scopic quantum tunnelling^ (MQT) have been observed 
in charge, flux, and phase qubits^. Renewed interest 
in MQT occurred after the recent discovery of MQT in 
high temperature layered superconductors^i^i^. The ob- 
served enhancement of MQT was attributed to the spa- 
tial structure of the tunnelling fluxon^i^. It is important 
to develop a theory of quantum electrodynamics in long 
(about 1 /im in stacks of J Js, and about tens of microns in 
low-Tc junctions) JJs where the spatial distribution of the 
gauge-invariant phase difference is crucial. In this prob- 
lem, the standard quantum mechanical approach (where 
the phase difference is associated with the coordinate of 
a quantum particle tunnelling through an effective po- 
tential barrier) becomes invalid, and a more advanced 
field-theoretical approach is needed^ii^. 

Here we consider a Josephson junction (JJ) driven by 
a dc current near its critical value and exposed to THz 
electromagnetic (EM) waves. In this configuration, as 
known for point-like contacts, the probability of MQT 
depends on the intensity and frequency u of the inci- 
dent EM waves. In contrast to the short-junction case, 
we predict several resonant enhancements of the MQT 
escape rate, when the frequency cu matches the eigenfre- 
quencies of the JPWs. We also propose a full quantum 
electrodynamical description of long JJs, to calculate the 
probabilities of absorption and emission of light by JPW 
quanta. 

In section II we derive the model and quantize the field 
of the gauge-invariant phase difference ip. In section III 
we consider the interaction of the quantized (p with pho- 
tons and calculate the transition rates of absorbtion and 



emission of light by JPW quanta. This allows us to find 
the mean values of occupation numbers of JPW quanta 
and the mean energy of the system, which is pumped 
by external THz radiation. In section IV we calculate 
the probability of photon-assisted macroscopic quantum 
tunnelling of the phase difference. 



II. SECOND QUANTIZATION OF THE PHASE 
DIFFERENCE FIELD 



A. Lagrangian formulation 

The geometry of the Josephson junction (JJ) under 
study is shown in Fig. [T] Two superconducting bars 
overlap a length D in the x direction. An insulating 
layer of thickness s, about several nanometers, is placed 
between these two bars. A supercurrent with density 
i flows through the junction in the z direction. The 
width L of the JJ in the y direction is of the order of, 
or less than, the Josephson penetration depth Aj, that 
is I — L/Xj < 1. The dynamics of the gauge- invariant 
phase difference ip(t, x, y) of such a junction is described 
by the action 
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In these equations, the x and y coordinates are normal- 
ized by Aj, the time t is normalized by l/tOp, where cUp is 
the Josephson plasma frequency, ic is the critical current 
density, and 
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where A is considered to be much larger than unity, A ^ 
1. The integration in £e['/5] is performed over a contour 
S around the junction's area, and the subscript n refers 
to the component normal to the contour E in the XY 
plane of the vector product of the magnetic field H and 
unit vector e^. 

The classical equation of motion for (/9 with the ac- 
tion ([1]) is the two dimensional Sine-Gordon equation 
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The surface term in the action ([T]) depicts the boundary 
conditions to this equation, 
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Representing the magnetic field in the form H = Hj + 
H"^, where H*^ is the external ac magnetic field and Hj 
is the field generated by the flowing current, we obtain 



dip 


l" 1 '''^' 


dx 


x=±dl2 2 ' 47ricAj 


dip 


^iil 


dy 


y=±l/2 47ricAj ' 


' = 


iD ^ D , L 
-^-, d= —, 1 = — 

IcA.l AJ Aj 



(4) 



where 



When H"^ = 0, the stationary solution to Eq. ([2]), 
corresponding to the lowest-energy minimum, does not 
depend on the y-coordinate. Below we assume that 
the ac magnetic field of the incident radiation is di- 
rected along the y-axis (see Fig. [T|). In this case, only 
the plasma waves in the x-direction are excited, and 
ip = ip{t, x). When D <C Aj, the field ip{t, x) only slightly 
depends on the s-coordinate, and the action ([T]) describes 
the dynamics of the particle in the washboard potential 
V{ip) — —cosip — jip, where j = i/ic- When j < I, 
this potential has an infinite number of minima, each one 
separated by a potential barrier of the order of \/l — p. 
The probability of quantum tunnelling from one mini- 
mum to the nearest minimum can be easily calculated in 
the semiclassical approximationii. When D > Xj, the 
spatial dependence of the field ip{t,x) is essential, and 
the problem of quantum tunnelling becomes more com- 
plicated. In the semiclassical approximation, the proba- 
bility of tunnelling can be written as^^ 



r = cjo\/ exp(-S) 



where B = 2S^/h, and iSe is the action, defined in 
Eq. II]), in imaginary time t = it calculated along clas- 
sical trajectories, and Wq is the oscillation frequency of 
the field ip{t^x) near one of the energy minima lpq{x). 
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FIG. 1: Schematic diagram of the Josephson junction. The 
wave vector of the externally applied polarized THz electro- 
magnetic wave is directed along the x axis, while its electric 
(magnetic) field is directed along the z (y) axis. 



In one of our previous papers^ we proposed an approach 
for calculating the tunnelling exponent B for a current 
/ = jd close to the critical value Ic{d) (which now non- 
linearly depends on d). Here we consider the effect of 
external electromagnetic radiation on the probability of 
tunnelling. 



B. Quantum regime 

We consider the interaction of ip with electromagnetic 
waves by perturbations. First, we quantize the field 
ip{t,x) near the energy minimmn ipoix) at H'^ — 0, find 
the energy spectrum, and then calculate the transition 
rates of the field (/?, from the ground state to its excited 
states and vice versa, due to the interaction with the 
electromagnetic field. The knowledge of the transition 
rates gives us the mean energy Eioj^P) of the field ip in 
the presence of an external radiation as a function of its 
power P and frequency to. Since the effective potential 
barrier decreases with the growth of E, the external ra- 
diation enhances the tunnelling. It is clear that a strong 
enhancement of the escape rate F should be at frequen- 
cies w close to the eigenfrequencies aj„ of the ip field. 
The tunnelling exponent B, as function of E, is found 
here using the approach described in Ref. |3- 

The static solution corresponding to an energy mini- 
mum satisfies the static Sine-Gordon equation 
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with the boundary conditions 
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The solution to this equation exists for currents / less 
than the critical value Ic{d)- If d ^ 4Aj, the current den- 
sity in the JJ is approximately constant and the function 
Ic{d) increases linearly with d; if d 3> Aj, the current 



flows near the junction edges and Ic{d) reaches the satu- 
ration value /™^'<^ — 4. In order to quantize ip we repre- 
sent it in the form 

(p{t,x) = ipo{x) +il;{t,x), 
where the operator ip satisfies the boundary conditions 



di; 
dx 



0, 



(7) 
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and expand the Lagrangian C in powers of ip. We intro- 
duce the momentum 



n(t, x) 



SC 



H^%^) 



dij{t,x) 
^hujpA — — — , (8) 



and require the standard simultaneous commutation re- 
lation 



ip(t,x),'!T{t,x') = ih6{x — x') . 



The Hamiltonian of the system, 7i, has a form 
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In these equations, the colons " : " mean "normal order- 
ing" and P is a differential operator of the form 



52 

V = -— -|-cos((/?o(x)). 
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In the interaction representation, the operators if) and tt 
can be written as 

^ = \/^E^(«^^*^l + e-^*^' (13) 
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where /z„, ipn are, respectively, the eigenvalues and or- 
thogonal eigenfunctions of the operator T), that is, 

.d/2 

Vlpn = fJ-n'ipn , / dx ■>pn{x)lpmix) = Snni , (14) 

J-d/2 

In Eq. ([T5| . 6|j and 6„ are the creation and annihilation 
operators of JPW quanta in the state n. Note that all 



/i„'s are positive when / < Ic{d) because the 1^90(2;) corre- 
sponds to an energy minimum. In terms of the operators 
6jj and 6„, the Hamiltonian Hq takes the form 



Ho = hujp 2J 



blL 



'^.n O'^fin 



(15) 



The Hamiltonian 7i' describes the self-interaction of the 
field ip. Since A :^ 1, Ti' can be considered as a per- 
turbation if the energy of the system (counting from the 
'vacuum' state corresponding to ipo{x)) E <C hUpA. In 
zeroth order, this energy is determined by the occupation 
numbers -/V„, and reads 



E ^ hujp^ ^/JI^ Nn . 



The correction to this result due to self interactions can 
be found via perturbation theory. 



III. INTERACTION WITH AN 

ELECTROMAGNETIC FIELD: ABSORPTION 

AND EMISSION TRANSITION RATES 

Now we consider the interaction of the field ip with 
electromagnetic waves, described by the vector potential 
A (we choose the gauge Aq = 0, divA — 0). Substituting 
ip — ipo + Ip into the action ([1]) and expanding it in a 
power series of "0, we derive the operator V describing 
the interaction of ip with the electromagnetic field^: 



V = 



cEj 



i-KicXjL 



dC^ [rotAxe,]„ . (16) 



Here we use the relation H^ = rotA/Aj because we mea- 
sure distances in units of Xj. The vector potential A in 
Eq. (|16p consists of two parts, describing both the incom- 
ing and outgoing radiation re-emitted by the J J. We as- 
sume that the incident electromagnetic radiation is fully 
polarized and propagates along the x axis, as shown in 
Fig. [TJ Below we measure the frequency co in units of Up 
and the wave length in units of Xj. In this case, A can 
be written as 
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where a^^ and atA are the creation and annihilation op- 
erators of a photon with wave vector k and polarization 
A, 

_ ojpXj 



is the ratio of the Swihart velocity ojpXj to the speed of 
Hght c, (jJk — |k|/w, and V is the volume of space (di- 
mensional) where the electromagnetic field exists. The 
first term in Eq. (|17p corresponds to incoming radiation 
(which is here considered as classical), where E^ is the 
electric field at frequency uj. The second term describes 
the photons appearing due to the interaction of the in- 
coming electromagnetic waves with the JJ. In this term, 
e"^ is the vector of polarization, which satisfies the equal- 
ity 



k-e^(k) = 0. 



(18) 



Substituting Eq. p?)) and the expansion (fT3|) for -0 into 
Eq. pop , and performing the surface integration, we de- 



rive: 
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where A^ is the London penetration depth, which is 
linked to ic and Xj by the relation (it is supposed here 
that s < Al) 



A^.= 



hc^ 
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h^(k) =kx e^(k), (22) 

and functions >f„(fc), Xn(fc), and S'(fc) are the following: 

>i^{k) = -I U„(d/2)e**'-'^/2 „ ^„( _ d/2)e-^'='^/2] , (23) 
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A. Spontaneous photon emission 

In first order of perturbation theory, there are only 
three possible processes: (i) spontaneous emission of a 
photon by the field if, (ii) induced photon absorption, 
and (iii) induced emissionii. Let us first consider the 
spontaneous emission, which is described by the operator 
Vq. In the initial state, we have the set of occupation 
numbers of JPW quanta, {N„i}, and zero photons, and, 
in the final state, one of these numbers, say iV„, decreases 
by 1 and one photon appears in the system. We neglect 
the effect of thermal radiation, proceeding to the limit 
T — > 0. The probability per unit time of such a process. 
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is proportional to Nn 



Following the standard approach^"^, we derive for the 
probability of emission of a photon with wave vector k 
and polarization A: 



dwl-\k)=N, 
where 
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(19) F,^(k) = S{ky) [kyh^{k)Xn{k.) + h^ik)K.n{k,)] . (27) 



We introduce spherical coordinates in momentum space. 
Performing the integration over k and the summation 
over A taking into account the relation p^ . finally, we 
derive 
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where m is a unit vector in momentum space. 

For not too long junctions, it is possible to obtain an 
analytical expression for ;/„. The wave length A of the 
electromagnetic radiation under consideration is about 

A -- c/up = Xj/v > Aj , 

since the typical value ofv ^ 3-10~^ <C 1. Therefore, for 
d -^ X, one can expand F,^(ci;„um) in Eq. (|^ in powers 
of V. Doing so, we derive in the lowest order 
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and h'2m ^ i^2m+i- The difference between z/„ with odd 
and even n comes from the symmetry properties of the 
JPW wave functions: ipni—x) = {—l)""tpn{x)- 



The value of 7„ gives us the radiation width of the nth 
level in units of cUp. When d ~ Xj, we have from ([50]) 
^^2™+! ~ v'^ and i^2m ^ v'^- Considering L, Aj ^^ 10^^ cm, 
Xl ~ 10~^ cm, and v ^ 3 ■ 10^^ we obtain 72m+i ~ 
lO-S-lO-"^ and 72m ~ lO^i^-lO^^''. Note that we do not 
consider here another possible mechanisms of dissipation, 
which can substantially increase the width of the JPW 
quanta energy levels. 



B. Induced photon absorbtion and emission 

Let us now consider processes of induced photon ab- 
sorption and emission. These two processes are de- 
termined by the operator V^xt- We denote by '"w„ 
("^uiJi ) the probability per unit time of creation (an- 
nihilation) of a quantum of the ip field in the nth state 
due to induced photon absorption (emission) . These two 
probabilities satisfy the following equality 



'w 



(-) 



Nn 



Thus, to first order in perturbation theory, the probabil- 
ity per unit time of induced photon absorption and also 
accounting for induced emission, m„ ' ='" Wn — '" Wn , 
does not depend on 7V„, and is only determined by the 
power and frequency of the external radiation. Making 
a similar calculation as for 
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We assume that the incident radiation has a Gaussian 
distribution with central frequency 6j and width 7, that 
is 
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where P is the radiation power per unit area and 
1 
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The probability Wn then becomes 
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where 
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In equilibrium, the probabilities Wn and Wn coin- 
cide. This gives rise to a relation for the mean values of 
the occupation numbers iV„: 
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The mean value E of the system energy (to zeroth order 
in 7i') then reads 



^ ^ A^PX^j ^ /„(0) 
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If the frequency band of the radiation source is large 
enough, that is, 7 3> 7n, we can easily perform an inte- 
gration in Eq. (|36p . As a result, the mean energy becomes 



UJp '^ 



>^«(wn'y)r 



r(LJ-u;„). (39) 



In the opposite case of near-monochromatic radiation, 
7 <C 7n, we should take into account that the energy 
levels of JPW quanta have finite width 7„ (in units of 
LOp). Replacing the delta function in Eq. ((36)) by 



(5(cj - w„) —t - 



In 
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and using T{lj — Cj) = 6{lu — Co), we obtain 
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where we take into account that 7„ ^ 1. Let us notice 
that formulas ([38]) . ((39|) . and (f40|) are valid only at not 
too high radiation power when E <^ hiUpA. 



C. Response of the junction to a wave packet 

Consider now the response of a JJ to a wide-band THz 
wave packet. We now assume that the central frequency 
of the incoming radiation uj is about tOp and that the 
width 7 of the wave packet is large enough. In this 
case, the first several energy levels of the system will be 
excited. The intensity U{u!) of light re-emission at fre- 
quency u! is given by the sum huip X^nA '^ndw^^^ (k), with 
dwl^ (k.) from Eq. ip^ . integrated over all directions of 
k. Taking into account the relation ((37)) for the mean val- 
ues of the occupation numbers N„, and replacing again 
the delta function in Eq. (|26p by a Lorentzian curve, we 
obtain 



f/(w) :^2PXlLJ2- 



{uj - UJnY + 72 



(41) 



The function U{llj), for relatively short {d — 2) and long 
JJs {d = 5), is shown in Fig. [2| The wave packet of 
the incident radiation has a central frequency a) = 1 and 
width 7 = 0.2 (in units of uip). In this case, the first two 
{d — 2) or three {d — 5) energy levels are excited. For 
short junctions, d < 1, the eigenfrequencies are, approx- 
imately 
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FIG. 2: (Color online) The frequency dependence of the in- 
tensity of re-emission, U{lj), calculated for d = 5 (red solid 
curve) and d — 2 (blue dashed curve) . Other parameters are: 
I/I^{d) = 0.98, Xl/\j = 4 ■ 10"^ and L/Xj = 2. The green 
dot-dashed curve corresponds to Gaussian distribution of the 
intensity of incoming radiation with central frequency Cj — 1 
and the width 7 — 0.2 (in units of tOp). 



into Eq. ((T|) and expanding the action in powers of ip, we 
obtain 



fro f rd/2 
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The last term in Eq. ([IS]) originates from the match- 
ing condition for the wave function $ of the quantum 
field i\) inside ($in oc exp(— 5^;/^)) and outside (<i>out « 
exp{—iEt/ti)) the barrier. The field tjj in Eq. ((45|) satisfies 
the equation 5B{E)/5%p = 0, that is 
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OT^ 2 6 
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with the following initial and boundary conditions 

= 0. 
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We seek a solution of the equation (j46| in the form 

V'lr.x) =^c„(r)Vn(a;). (47) 



For increasing values of d, wi tends to wq, and when 
d > 4, we have 



UJQ 



t^i «C cj„ , n > 1 . 



(43) 



The relation uji w loq is essential for the properties of 
macroscopic quantum tunnelling in JJs. Namely, in this 
case we have two channels of tunnelling, corresponding to 
fluxons arising near the junction's edges. This situation 
is considered in the next section. 



IV. PROBABILITY OF PHOTON-ASSISTED 
TUNNELLING 



Now we calculate the probability per unit time of quan- 
tum tunnelling to another vacuum state, stimulated by 
external electromagnetic radiation, using the approach^ 
proposed in one of our previous papers. In the semiclas- 
sical approximation, we can consider the quantum field 
■0 as a classical field ip(T,x) in imaginary time t — it. 
The probability, F, then readai^ 



(44) 
where Se is the action ([T]) in imaginary time. Substitut- 
ing 



Multiplying Eq. (I46p by ipn and performing space inte- 
gration and using Eq. (fT4|) , we obtain the system of equa- 
tions for Cnir) 
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with initial conditions 

c„(0) = c„(to) = 0. (49) 

Here, the dot means "imaginary time derivative" , and 
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The tunnelling exponent B{E), Eq. ([45|) . can be ex- 
pressed as 



B{E) = A / dr 



1 V U^^^ 
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When the current / is close to the critical value Ic{d), 
we have /j,o ^ 1 and c„ ^ 1. So, we can neglect all 
terms in the right-hand-side of Eq. ([48|) . except the first 
one. Our analysis shows that when d > 4, //i « /iq, and 
we have the following relation for the eigenvalues of the 
operator V 



(p{T,x) = ipo{x) +Tp{T,x) 



Mo ~ Ml < Mn : n> 1, 




FIG. 3: (Color online) The potential V{ao,ai) calculated 
when d = 6, I/Ic{d) = 0.98 {X/uoi = 1.18). A particle tun- 
nels from its initial position near a^ — 0. The Qi are collective 
coordinates for the tunnelling fluxon. The curves correspond 
to three possible imaginary-time trajectories of the particle: 
al {rj) and a^ {r)) (see text below). Note, that the real-time 
potential equals to — V(ao,ai). 



In this case, co,Ci <C c„ (n > 1), and we can 
only the first two equations of the system ([48 
c„ = for all n > 1 (for details, see Ref. [3)- 
We now introduce new variables 



a.(.) = ^^^f^, ^ = 0,1, 



where 



V = y/lM)T , Ui = U^. 



The system of equations 
dPao 



(Pai 
drf' 



ao 



(3) 
Oil ■ 

takes the form 
-- (a^ + a-^) , 



consider 
I, taking 



(52) 



(53) 



(54) 



Aai = — 3uoi aotti 



where 



Mo 



■, "01 






The system ((M)) has the first integral 



%4 



dao\ 1 fdai\ 

drj J uoi \ drj J 



where we introduce a potential 



T^(ao, ai) = a'l + 2>aQa\ — an of 

uoi 



(55) 



E 
(56) 

(57) 



to 
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FIG. 4: (Color online) The frequency dependence of 
b{Q)/b{e{uLi,P)) — 1, calculated for d — 2 (upper graph) and 
d — 5 (lower graph). Other parameters are I/Ic{d) = 0.98, 
L = 2Xj, and v = 1/30 for both cases. Solid curves 
correspond to monochromatic incoming radiation, while the 
dashed curves describe the response on wide-band THz ra- 
diation with 7 = 5- 10~^ (in units of ujp). The radiation 
power P is the same for all curves and is chosen such that 
max(e(a;, P)) — 0.3eo ~ 2/45 for wide-band radiation. 



we have. 



Taking into account the initial conditions (j49 
at the turning points: 



V{ao{^T), ai{y/JI^T))\^^^^_^^= 



where 



e = 






< £ < £0 = ^ 



(58) 



(59) 



The equation ((58)) defines the value of tq as a function of 
system energy E. 

Thus, we reduce the problem of quantum tunnelling 
of the field (p to the problem of tunnelling a quantum 
particle in two dimensions, where the a.i's play the role 
of the particle generalized "coordinates" . The potential 
V{ao, ai) is shown in Fig. [31 

When d > dc{I,e) « 4, there are three solutions of 

the system of Eqs. ([54)) with the conditions ()58p . al (77) 

and a\ (tj), which are characterized by the following 
relations 



,.(0) 



(r,)=0, a[ \v) 



,.(+) 



iv)- 



The trajectories a- (rj) and a\ (rj) are shown in Fig. [3] 



S+), 



„(-) 



The solution a\ (rj) {a\ (?])) corresponds to the forma- 
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tion of vortex (antivortex) nucleus at left (right) junc- 
tion's edge, while the solution a,- (ry) describes the tun- 
nelling of (p as a wholeS. The minimum of B{E) corre- 
sponds to the solutions al (77) . The tunnelling exponent 
then reads 



B{E) = 



5wg 



bis) ., 



(60) 



where 



bis) 



15 
16, 



rvo 

/ dr] [ag -I- Saotti — 2 e] , 
Jo 



Vo = VA^oTo ■ 



(61) 

Note, that we should multiply the probability ^{E), 
Eq. (j44|) . by a factor of 2, since we have two channels 
for tunnelling. 

When d < dc{I,£), all three solutions coincide, 

Q^i iv) =^ Q!i iv) — ^1 iv) — Oj Eind the second equa- 
tion of the system ([5^ becomes trivial, while the first 
one can be easily integrated. As a result, we obtain 



m 



15 
16 



2(e) 

da- 



2s 



Qi(e) 



V«'(l-«) 



5(0) 



(62) 



where 01^2(2) are the smaller and larger positive roots of 
the cubic equation 



a^{l-a) 



0. 



The analysis of the tunnelling exponent B on the junc- 
tion's width d and current / was carried out in one^ of 
our previous papers. Now we are interested in the ef- 
fect of electromagnetic radiation on B{E). Using formu- 
las (gni), dSH), dSJ), and dSl]), we calculate the depen- 
dence of 6(0)/6(e((I',P)) — 1 as a function of radiation's 
central frequency (D, for short d < d^, and long d > dc 
junctions. The results of the calculations, both for wide- 
band (7 ^ 7n) and monochromatic (7 <C 7n) radiation, 
are shown in Fig. [H It is clear that we have several res- 
onances at frequencies Cu — Un — y/Jhi (in units of Up). 
When 7 = 0, the resonance peaks are very narrow, and 
one can switch the JJ to the resistive state (B = 0) at 
small radiation power, ii uj — LUn- Note, that the con- 
dition Lj — LUn can be achieved by changing the applied 
dc current /, since ujn depend on /. In other words, if 
the frequency of the incoming radiation lies near one of 
the LUn , one can observe a resonance behavior of the tun- 
nelling exponent i? as a function of the dc current. 



For relatively short junctions, d < 1, we have several 
well separated peaks even for wide-band THz radiation, 
as it can be seen from Fig. 21 while for d > dc, we have 
luq ~ LUi <^ LUn, and the first two peaks can merge into 
a single peak. Note, that the inequality cjq ~ wi <C LUn 
is valid for not too high junction's width d < 20. In 
the opposite case we should consider the large number of 
equations in the system (|48p to calculate the tunnelling 
exponent. 

V. CONCLUSION 

In conclusion, we proposed a quantum field theory for 
Josephson plasma waves interacting with external elec- 
tromagnetic waves. We also calculate the macroscopic 
quantum tunnelling of a fluxon, stimulated by THz light, 
in a long Josephson junction driven by a dc current. The 
probability of absorption and emission of light depends 
on the current and the length of the Josephson junction. 
The MQT escape rate shows several resonance maxima 
as a function of the frequency, corresponding to eigenfre- 
quencies of Josephson plasma wave quanta. This could be 
potentially useful for a variety of superconducting quan- 
tum THz devices. Classical THz devices are discussed in 
Ref. Hi. 
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